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This study analyzes thermoelectric properties of a one-dimensional random conductor which shows 
localization effects and simultaneously includes resonant scatterers yielding sharp conductance res¬ 
onances. These sharp features give rise to a distinct behavior of the Seebeck coefficient in finite sys¬ 
tems and incorporate the degree of localization as a means to enhance thermoelectric performance, 
in principle. The model for non-interacting electrons is discussed within the Landauer-Biittiker 
formalism such that analytical treatment is possible for a wide range of properties, if a special aver¬ 
aging scheme is applied. The approximations in the averaging procedure are tested with numerical 
evaluations showing good qualitative agreement, with some limited quantitative disagreement. The 
validity of low-temperature Mott’s formula is determined and a good approximation is developed 
for the intermediate temperature range. In both regimes the intricate interplay between Anderson 
localization due to disorder and conductance resonances of the disorder potential is analyzed. 


I. INTRODUCTION 

In the context of energy harvesting, research on ther- 
moelectrics has acquired renewed momentum, as efficient 
thermoelectric materials may provide a means to con¬ 
vert heat into electrical energy at relatively low main¬ 
tenance and reliably as no engines with moving parts 
are involvecP®^. Thermoelectricity constitutes a stan¬ 
dard part of transport theory: electricity and heat are 
connected, for instance, by the Seebeck coefficient S 
which relates a temperature difference AT with a volt¬ 
age difference AU as AV = SAT assuming open circuit 
condition^^H^. An important quantity for applications 
is the figure of merit Z, a measure for the efficiency of 
the energy conversion for a material acting as a ther¬ 
moelectric device. The dimensionless parameter ZT is 
defined as ZT = aeiS'^T/n, where aei (k) denote the 
electrical (heat) conductivity. Much effort is devoted to 
enhanciM ZT for which values beyond 1 are only rarely 
reportecF^*^. Strategies for improvements have turned to 
nanostructuring of materials which reduces the phonon 
heat conductivity and increases a and S by adjusting in¬ 
ternal properties as mo bile charge carriers confined in a 
narrow enemy rangd^*2H121 Alternatively, also correlated 
systemalSHlSl and low dimensionality have been consid¬ 
ered as a way to strongly suppress the phonon heat con¬ 
ductivity and enha nce th e electrical conductivity as well 
as the thermop owei * ^^ 

In the low-temperature limit the Seebeck coefficient is 
given by Mott’s formula, 

c_ 7r2fc|r dlog aei 

which relates S to the energy dependence of the electri¬ 
cal conductivity a^i at the chemical potentiaP^^. This 
formula indicates that S is also a measure of the energy 
dependence of the conductivity. This aspect has been, 
for example, been emphasized by Mahan and Sofo who 
suggested to optimize the thermoelectric performance by 


using devices where the conductance has delta-peak like 
structures in the vicinity of the chemical potentiaP^. 
Such kind of structures are naturally obtained in systems 
with conductance resonances. 

In our study we focus on disordered one-dimensional 
systems with disorder where we investigate the effect of 
Anderson localization and, through special design of our 
model, also the situation of mobile charge carriers near 
localization. Such a system can be realized by a model of 
randomly positioned impurities, barriers of given speci¬ 
fications. In general, all carriers are localized in such a 
system. However, if these scatterers are identical, they 
can develop conductance resonances at specific values of 
energy. In combination with the Anderson localization 
effect such resonances show interesting features which, 
in principle, could be used to design improved thermo¬ 
electric devices. Note that Anderson localization in gen¬ 
eral has been investigated by many groups for a variety 
of reasoiP^lt^. Using the transfer matrix formulation by 
Landauer and Biittiker we will analyze different regimes 
of the system analytically as well as numerically. This is 
possible if we ignore the interaction among the electrons. 
While transport properties of electronic states near the 
mobility edge of systems displaying Anderson localiza¬ 
tion physics are often effectively modeled, for instance 
by a variable range hopping model, we have not to resort 
here to any effective model. Rather we will benefit from 
the fact that the transfer matrix approach enables us to 
deal with various aspects of a disordered one-dimensional 
system analytically. While we are mainly interested in 
the basic behavior of thermopower in our special type of 
random one-dimensional model, we will also briefly ad¬ 
dress statistical features of a finite random system for 
comparison with related discussion^. 

II. MODEL 

We first introduce a model of a disordered one¬ 
dimensional system of non-interacting electrons. Con- 
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figuring the model as a random array of well localized 
elastic scatterers will allow us to deal with the effect of 
Anderson localization to a large extent analytically by 
using the transfer matrix approach. 


A. Transfer Matrix Method 

The Landauer-Biittiker formalism of transfer matrices 
is undoubtedly the most successful method to discuss 
transport properties of one-dimensional mesoscopic sys¬ 
tems, if only elastic scattering is involved which conserves 
the electron energji^SlMI Inelastic scattering processes, 
such as electron-electron or electron-phonon scattering, 
are omitted for simplicity. 

The effect of scattering of an electron at a potential 
limited to certain range in space is encoded in the so- 
called transfer matrix T, 



This matrix relates the in and out-coming wave functions 
at the left hand side, A and B, with the out and incom¬ 
ing wave functions at the right hand side of the given 
spatial range of the potential, C and D, see Fig.[^ in the 
following way 



( 3 ) 
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FIG. 1: Notation of the in- and out-going wave 
functions for a given spatial range of the potential 
described by T. 

The system resistance introduced by backscattering of 
electrons is given by 


IrP 

Rsys = 1 (4) 

where i?o = /i/2e^ and 1 = |rp-|-|tp. The measured resis¬ 
tance includes also the contact resistance Rc = h/2e^ = 
i?o which gives then 


R = R,y,+R,= ^. (5) 

The full conductance is the inverse of i?, 

G = Go\tW ( 6 ) 

where Gq = 2e^//i is the conductance quantum. 


We model our system by well localized potential barri¬ 
ers. For mathematical simplicity we assume box-shaped 
potentials as sketched in Fig. which can be easily 
parametrized by the height Vi and the width 5i and the 
distance h to the potential on the left hand side. 



Iq li Iji 

FIG. 2: Sketch of the potential landscape with its 
underlying parameters. 


For a system with N barriers with a given set of pa¬ 
rameters {{li,Vi,6i)}i=i^,,,N we can compute the corre¬ 
sponding transfer matrix T by splitting it into elemen¬ 
tary transfer matrices, Tp for the propagation between 
the barriers (scattering free) and T) for the scattering 
part of a potential barrier, 

f = fp{lN) ■ f,iVN, Sn) ■ fp{lN-l) ... RiVi, Si) ■ fpilo) , 

( 7 ) 

which corresponds to a simple sequential product of 
transfer matrices. The propagation transfer matrix spec¬ 
ifies only the change of the phase between barriers and is 
given by 


fp{l) 


/^ikl Q \ 

0 ’ 


( 8 ) 


whereas the impurity transfer matrix is more complex. 
It can be computed by solving the one dimensional 
Schrodinger equation and imposing the continuity condi¬ 
tions at the interfaces. For our box potentials this yields 
after some straightforward calculation. 


f(VS\= (+ f sm{qS) f sm{qS) \ 

^ -fsin((z7 COs{qS)-fsm{q6))^ 

( 9 ) 

where k = \/2mE/h, q = ^2m{E — V)/h and 


e = 




( 10 ) 


By combining Eq. dWI ) we can compute G for any 
configuration of barriers in a system. The link to ther¬ 
moelectricity is obtained through the Seebeck coefficient 
S given by the Gutler-Mott formula. 


SG{T,fi) 


1 JdE {E-f,)GiE) (-§1) 
JdEG(E)(-§l) 


( 11 ) 


where / denotes the Fermi-Dirac distribution function. 

This expression for the thermopower includes the con¬ 
tact resistance and differs from the thermopower of the 
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wire only. However, it has been showiP^ that this dis¬ 
crepancy vanishes in the limit of infinite scattering cen¬ 
ters. Therefore, we do not specifically distinguish be¬ 
tween these two kinds of Seebeck coefficients and use the 
formula in Eq. (11). 


It turns out that the former is more easily accessible in an 
analytic approach and we will for the time being focus on 
this approach. However, later we will compare the two 
averages in a numerical discussion. 


B. Averaging Procedure 

The Seebeck coefficient So depends on the configura¬ 
tion {(Zi, Vi, Si)}i, which constitutes a too large number of 
parameters. Therefore we turn here to an averaging over 
many configurations, assuming a certain self-averaging 
for large enough systems. The averaged Seebeck coeffi¬ 
cient (S) is defined by 

iS)=J (^dkdV,ds}j P{{ik,V,S,)},)SG , ( 12 ) 

with the probabilistic parameter distribution function P. 
This definition is very general and a direct analytical 
evaluation of (S') is challenging. 

We will work with an alternative form of averaged ther¬ 
mopower. The basic idea is to average first the conduc¬ 
tance G instead of S, see Fig.[^ and secondly, to calculate 
the thermopower by applying Eq. 0 on the average of 
G. Later we will discuss this scheme by comparing the 
two averaging procedure using numerics. 



FIG. 3: Schematic diagram of the averaging process. 

We distinguish here two averages the harmonic and the 
arithmetic. The harmonic average of the conductance, 
{G)h, is given by 



while the arithmetic one, {G)a, has the form, 



C. Distribution of scattering barriers 

The energy dependence of the conductance is an inter¬ 
ference effect. Whenever the width of a barrier or the 
distance between two potential barriers matches a multi¬ 
ple of the electrons wave length, constructive interference 
and perfect transmission through this part of the system 
occurs. Due to the random distribution of the scatter- 
ers such resonances are washed out. Here, however, we 
would like to introduce a certain distribution which al¬ 
lows to see perfect transmission. For this purpose we 
assume that the different parameters are distributed in¬ 
dependently. Thus, we can assume the product form, 

P{{ih, Si, Vi)}i) = Pi{5j, Vj) Pp{lj) 

3 

= X{m 3 ) PsiSj) Pv{V,) , (15) 

3 

where 

Piili) = const , (16) 

Pvm=6iVo-Vi) , (17) 

using the (5-function (5(...) in Py (not to be confused 
with the impurity width S). Inter-barrier resonances are 
completely wiped out by this distribution and the barrier 
height of all scatterers is fixed. This will lead generally to 
a localization and a vanishing conductance for infinitely 
large system. However, our model has built-in the possi¬ 
bility for a recovery of the conductance by perfect trans¬ 
mission, if all barriers satisfy the resonance condition si¬ 
multaneously. This can be reached for a — 0 where we 
find conductance peaks for specific resonance energies of 
the electrons. These perfect conductance resonances are, 
however, reduced and washed out with growing standard 
deviation a. 

Thus, perfect transmission can now only be reached 
within the barriers and be wielded using the standard 
deviation a of the widths Si. 

III. ANALYSIS OF THE SEEBECK 
COEFFICIENT 

The above approximations allow us to continue our an¬ 
alytical study of thermoelectricity and will give us good 
insight in the basic properties in various regimes and the 
important parameters. 
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A. Computation of S{G) 
1. Averaged Conductance (G) 


We start with the averaging of the conductance and 
use, as announced above, the harmonic average which 
gives us a simple analytical result. The harmonic average 
of G defined in Eq. (131 is given by 


(G) 


2Go 

1-f 


11 


(19) 


both, the distribution functions Pp^i and the transfer ma¬ 
trices Tp i. The barrier heights Vi are already fix to be 
Vo according to the distribution function Py. The pa¬ 
rameter N is the number of barriers (impurities) in the 
system. For the explicit derivation of this re sult as well 
as the definition of Fp^i we refer to Appendix 


A1 


The explicit form of the conductance reads, see Ap¬ 
pendix |A2[ 


where Fp and Fi are 3x3 matrices which incorporates 


(GKx) = 


2Go 


2Gn 


+ (l + + (1 - “‘t'gA) ) 


N 


1-I-|fi(a;, A,cr)| 2 JV 


( 20 ) 


with the dimensionless energy x = FjVo-, the barrier pa¬ 
rameter A = 2mVo5 ‘^and |tip the mean transmission 
probability for a single barrier. 

The behavior of (G) is displayed in Fig. [41(a ) for various 
values of a/5. As is obvious from Eq. ( |20| ) resonances 
appear at energy values Xn = l-l-(n7r)^/A with the integer 
u > 1, yielding perfect conductance for ct = 0 (absolutely 
identical barriers, randomly located), i.e. {G)(xn) = Go- 
With increasing a the conductance shrinks. Sufficiently 
far from the resonance energies the conductance vanishes, 
manifesting Anderson localization. Note that in these 
intermediate energy ranges the conductance can be safely 
approximated by 

-1 

= 2GoMx,\,a)\^^ ( 21 ) 

11 

for ^ 1, which allows us to define a localization length 

(G) = 2Goe-^/^ , (22) 

where L is the system length with Lui = N {ui is the 
constant impurity density) and £ = —{2ni In |ti|)“^. 

The conductance resonances allow us now to approach 
an energy regime where the Anderson localization length 
can become comparable to the system size and the carri¬ 
ers become delocalized. Additionally, the regions around 
the resonances are interesting for thermoelectricity due 
to the strong energy dependence of the conductance, as 
the Mott formula suggests. 


(G) « 2Go 


(4c) 


N 


2. The Thermopower S^g) 

In Fig. |^b,c) we show the Seebeck coefficient S(^c) 
a function of temperature in two regimes which display 


characteristic behaviors. The data are given in units of 
'S'o = fcs/e « 86.2 [TV/K and To = Vb/^B- 


We use A = 9 and N = 1000 as the system parameters 
and check the temperature range 10“"^ < T/To < 10“^. 


In Fig. j^b) we set the chemical potential slightly be¬ 
low the lowest conductance peak which yields a negative 
Seebeck coefficient. Under this condition the resonance 
closest to the chemical potential dominates the behavior 
of the thermoelectric effect. Using a logarithmic tem¬ 
perature scale we see a pronounced peak in S(^g)- Its 
magnitude and the location are almost unaffected by a. 
Since the integrals in Eq. (II) are dominated by a sin¬ 
gle resonance, the reduction of the denominator for non¬ 
vanishing a is compensated by the analog reduction of 
the nominator. Only for a/5 > 0.2, the Seebeck peak 
starts to diminish and another peak arises at higher tem¬ 
perature which originates from the next conductance res¬ 
onance. 


In Fig.j^c), the chemical potential is between the first 
and second resonance but closer to the lower than the 
upper one. As a consequence, at low temperatures Si^q) 
is first dominated by the lower resonance which yields 
a positive thermopower. With growing T the influence 
of the broader upper resonance appears and eventually 
turns even negative. Interestingly, in this higher 

temperature regime even a rather pronounced peaks can 
appear for growing a. The reason lies in the different 
suppression of the resonance peaks in (G). The upper 
resonance is more slowly reduced, see Fig. j^a). 
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T/To 



FIG. 4: (a) The harmonically averaged conductance for 
ajb = 0,0.03,0.06 (solid, dashed and dotted line) at 
N = 1000 and A = 9. Based on this conductance we 
computed the thermopower S(g) for (b) 

P-i = = 1-6 (which is below the resonance energy 

El = 1 + 7r^/9 « 2.1) as well as (c) {I 2 = 2.3 (which is 
above the resonance energy Ei) where Sq = ks/e and 
To = Vn/kB. 


B. Analytical aspects of the thermopower S(G) 


which is valid at low temperatures where the ther¬ 
mopower is dominated by low energy excitations of the 
charge carriers at the chemical potential /r. We use now 
Eq. (21) and obtain with Eq. (23), 


^(G> = - 




3 e 


dE 


= -s< 


2Tr^ NkeT 
'IT Er 


where we defined a characteristic energy scale Ec by 


(24) 


d\og\ti 


dE 


d\ti\ 1 


dE \ti 


1 


(25) 


Note that E^. contains the information about the conduc¬ 
tance (G) as well as the position of the chemical poten¬ 
tial. In Eq. (24) we observe the surprising result that 


in the Mott regime the Seebeck coefficient is propor¬ 
tional to N, a result which had been previously obtained 
numerica lljE3. Thus, the higher the number of impurities 
the better the thermoelectric performance is found. 

The question arises up to which temperature the Mott 
regime is a valid approximation. Therefore we define a 
critical temperature Tcrit which we consider as an up¬ 
per bound for the Mott regime. We calculate S^c) with 
Eq. (11) as well as Mott’s formula and determine Tcrit, 


the temperature at which the relative discrepancy ex¬ 
ceeds 5%, as a function of N. The result is displayed 
in Fig. as a double-logarithmic plot, showing that 
Tcrit OC 



FIG. 5: Numerical calculation of the critical 
temperature as a function of the number of 
impurities N (dots) and a 1/iV fitting function (solid 
line) for p, = fi/Vo = 4,N = 1000, A = 9. 


1. Validity of Mott’s formula 


The standard approximation for the thermopowe 
Eq. (11) in the low-temperature regime is Mott’s formul 


We can argue on this A-dependence of Tcrit by the 
following discussion. Since in the integrals only the en¬ 
ergy range close to /i is important at low-temperature, 
we may expand ti{E) around yt assuming a weak energy 
dependence. 


%>=- 


7r2fc|T 91og(G) 


(23) 


ti{E) « fi(^) 


1-h 


E — yt 

Ec 


3 e 


dE 


(26) 
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which yields 


This inserted into Eq. (Ill gives 

^0 Ki 


%>=- 


ksT Ko ’ 


with 


Kr.=J dEiE-fir{G){E) 
{keT) 


I P + OO 

-A / dx x" 

4 ./__ '■'■"^2 


2x'y 


cosh (x/ 2 ) 


(27) 


(28) 


(29) 


where 7 = NksT/E^ and A = 2Goe^^*"l*i^^^L Obvi¬ 
ously, the low-temperature limit is only well defined if 


7 = 


NknT 1 
^ 2 ’ 


Er_ 


(30) 


such that the upper limit can be estimated as 


where AG = E^ — p. is the distance between the reso¬ 
nance energy and the chemical potential, and A « 0.6324 
(see Appendix B 1|B 2 for a detailed derivation). From 
this expression it becomes obvious that the sign of S(^c) 
depends on the sign of AE, i.e. whether the chemical po¬ 
tentials is above or below the resonance. As Fig.j^shows, 
the (5E-expansion matches perfectly the high tempera¬ 
ture limit of 5 ^( 3 ) (T) and is also able to include qual¬ 
itatively the thermopower extremum which is invisible 
within Mott’s formula. The extremum can easily be dis¬ 
cussed by using the following approximation, 


S(g)(T)^-S, 


( \^E\ 

[ksT 


A 5E^ 

8 (/cbT )2 


■sign(AF;). (34) 


The computation of the extremal point is then straight¬ 
forward. 


ksEmax — . 


A 5E^ 


(35) 


4 lAGl’ 


^ keN ’ 

as obtained numerically. It also indicates that the valid¬ 
ity of Mott’s formula is restricted to the very low tem¬ 
perature regime which shrinks with the system size, if 
localization plays a role. In order too discuss the high 
temperature regime and the influence of the conductance 
resonances we will have to go beyond Mott’s approach. 


Note that S^ax cc A, indicating that localization tends 
to increase the thermoelectric performance due to the 
sharpening of the resonance peaks. It is interesting to 
remark that within our approach the product 


\Emax ■s„ 


|ae| 

2 e 


(37) 


becomes independent of N, the number of impurities. 


2. Beyond Mott’s formula: the 5E-expansion 


Mott’s formula is inappropriate to treat the extrema 
of thermopower observed in our calculation displayed in 
Fig.|4l:b,c ). We resort therefore to a different approach. 
For simplicity, we focus on a single conductance peak 
at the resonance energy for a = 0 and expand the 
expression in Eq. (20 1 near E„, 


(G) 


2Gq 


i + (i + r„(E-F;„)2) 


N 


(32) 



which has a width at half maximum of SE « -^/s/lVN. 
Note that SE oc indicates that localization yields 

sharper resonances as states beyond the resonance energy 
tend to be localized. 

Since SE —?> 0 as A —>■ 00 we will use SE as an expan¬ 
sion parameter and express the Seebeck coefficient as 


FIG. 6 : The low temperature expansion (Mott) and the 
SE expansion in direct comparison with numerical data 
of S^G) for A = 9, A = 1000, h/Vq = 1.7. 


*5(0 


1 

-^0 


■ AE A SE^ 

k^ ~ 8 [keTY 


AE \ 
2kBT) 


,(33) 


3. Several conductance peaks 

The statements based on the JA-expansion can be ex¬ 
tended to the multi-peak scenario. Let us focus on low 
temperatures, ksT < A„+i—A„, where the thermopower 
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FIG. 7: (Color online) (a) We plotted the thermopower S multiplied by such that the extrema are equally 
pronounced for A = 9, = 8000. In (b) we plotted the thermopower for the temperature T/Tq = 0.1 which 

demonstrates that the SE expansion is in good agreement with the numerical data as long as T > Tcru which is true 

for the lowest two conductance peaks. 


shows a single extremum. It corresponds to the conduc¬ 
tance peak which dominates the integrals in Eq. (11) and 
lies closest to the chemical potential fj,. 

In Fig. I^a) we plot the numerical evaluation of the 
Seebeck coefficient S' as a function of T and /i. For illus¬ 
tration purpose we multiply S with the factor such 
that the extrema are equally pronounced. The green line 
describes the temperature T^ax of the extrema as a func¬ 
tion of /i, according to Eq. (35). The locations of diverg¬ 
ing Tmax mark the resonance energies. In addition, we 
display a cut along T/Tq = 0.1 to show S as a function of 
/i in Fig. IZtb)- We can here compare the numerical and 
the approximate result of Eq. (36) using always the clos¬ 
est conductance peak. The numerical data support the 
validity of a single peak picture, if the chemical potential 
is in the vicinity of a conductance peak. Then the ther¬ 
mopower shows a linear behavior with a sign change as 
a function of in agreement with the analytical results 
for T > Tmax (for T < T^ax tbe JE-expansions breaks 
down). Roughly in the middle between two resonance 
energies, there is an abrupt change of sign of the ther¬ 
mopower which corresponds to the switch in the domi¬ 
nance of the adjacent resonances. 

The temperatures T^ax of the individual conductance 
peaks are qualitatively in agreement with the numer¬ 
ical data, although there are discrepancies due to the 
influence of neighboring resonances, which is more pro¬ 
nounced for chemical potentials above the dominant res¬ 
onance. 


C. 5i?-expansion of the ZT-value 

The figure of merit ZT of a thermoelectric device de¬ 
scribes its efficiency. It is given by 

ZT = — (38) 

^el T ^ph 


with the electrical conductivity aei, the electronic and 
phonon heat conductivity Ke and Uph, respectively. 

If we simply neglect Kph, the JE-expansion of ZT can 
be straightforwardly evaluated (see Appendix |B 3[ ) which 
in leading order in SE gives, 


ZT 


8 

A 



(39) 


This result is in agreement with numerical calculations 
of ZT as long as T > Tmax- The apparently unlimited 
range of ZT values is a consequence of Kpt = 0, be¬ 
cause the dependence on SE~^ originates from Kei oc SE^. 
By taking the phonon heat conductance into account the 
ZT-value stays limited for all SE and AE. We parame¬ 
terize the phonon heat conductance as 

«p/.(r)=(^) ^Tp,,{T) = aTpniT) , (40) 

where a = 3.356 • 10“^"* W/K and TphiT) describes the 
temperature dependence. A straightforward calculation 
leads to the maximum of the ZT-vahxe with respect to 
SE and AE for given temperature T, 

/ T \ 

(ZT)max = niax ZT « 7 ( -^^ ) , (41) 

se,ae ' \Tph{T) J ’ ^ ’ 

where 7 = 7.696. The optimal value for SE and AE 
which maximizes the ZT-vahie are given by 


AE s 

s ±3.24 

■ ksT , 

(42) 


s 0.417- 


(43) 


This optimal value for AE is qualitatively in agreement 
with the results by Mahan and SofcP^ who performed 
similar computations for systems described by the Boltz¬ 
mann transport theory. They found an optimal value of 
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FIG. 8 : A comparison of the arithmetic, harmonic, renormalized averaged conductance and the numerically 

calculated conductance. 


AE = 2.4A:bT for a 5-function shaped transport distribu¬ 
tion function (mobile density of states), i.e. for 6E = 0. 
This discrepancy arises from the upper boundary Go of 
the conductance G which prevents a 5-function shaped 
conductance G and yields a finite optimal width of the 
peak, 6E ^ 0, and a larger coefficient for AE. 


IV. NUMERICAL SIMULATIONS 

The Seebeck coefficient derived via the harmoni¬ 
cally averaged conductance {G)h can be treated analyt¬ 
ically in various limits, in particular, also using the SE- 
expansion. However, the harmonic average of the con¬ 
ductance {G)h takes only the broad resonance peaks due 
to resonant scattering in the impurity barrier potential 
into account. It underestimates the contributions from 
the narrow peaks due to resonant scattering between the 
impurities, as we illustrate in Fig.|^ We compare the dif¬ 
ferent averages with the numerical calculation of the con¬ 
ductance for a concrete configuration which involves nu¬ 
merous narrow resonances. The arithmetic average does 
much better in this respect. 


A. Validity of the averages for S 

In order to understand this discrepancy qualitatively 
we define the renormalized averaged conductance (G)r. 
This averaged conductance is based on the harmonic av¬ 
erage {G)h, but includes the influence of the very narrow 
peaks between the broad resonances. We introduce the 
mean distance A between neighboring narrow peaks and 
their mean width w. Since the harmonic average is dom¬ 
inated by the smallest values, {G)h can be regarded as 
the minimal conductance, consistently with Fig. To 
this minimal value we should add the contributions of 


the narrow resonances, each one with the mean weight 
w = G ■ u! ■ (Go — {G)h) where G ~ 1 is a numerical 
factor. The renormalized harmonic average {G)r takes 
all those small weights into account by distributing them 
uniformly, 

(G), = (G)„ + g|(Go-(G)^). (44) 

Indeed Fig. shows that the renormalized average {G)r 
is much closer to the arithmetical average (G)q than the 
harmonic average. Thus we would expect that this also 
leads to numerically more accurate results for the Seebeck 
coefficient. 

For this discussion let us consider the numerically de¬ 
termined Seebeck coefficient {S) and compare it with our 
previous averaging procedure based on the harmonic av¬ 
erage of the conductance. We analyze the same regimes 
as in Sect |III at A = 9, V = 1000. We start with 
= 1-6 which yields a single thermopower extremum 
in the chosen temperature range. For the numerical av¬ 
erage we used 50 samples randomly picked according to 
the distribution function in Eq. (15), calculated the con¬ 


ductance and the Seebeck coefficient by means of the 
Cutler-Mott formula Eq. <0 which then is averaged over 
all samples. We can compare now the results for (S) in 
Eig. I^a) with the analogous results for S(^g) in Fig. 0i>) 
obtained with {G)h- Obviously, the qualitative features 
are identical: there is a maximum at roughly the same 
temperature and also the trends as a function of the stan¬ 
dard deviation a agree. The largest difference lies in the 
absolute magnitude, S(^g) is roughly three times the size 
of (5). 

The discrepancy can be accounted for by analyzing 
the consequences of the renormalized average (G)^. A 
straightforward calculation leads to a renormalized width 
of the conductance peak, 


5K — 

UJ-Jren — 


SE 


\/l - 2Gf ’ 


( 45 ) 
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FIG. 9: Numeric results of (S) for (a) fi = h/Vq = 1.6 
(slightly below the hrst resonance) and (b) jl = 2.3 
(between the hrst and the second resonance) with 
N = 1000 and A = 9, see Fig. 


which can used to determine the inhuence on the ex¬ 
tremum by using Eq. ( |35lMt , 


T 

rpren max 

max ‘2.C ^ 


qren _ 

'^max 


(l - 2C^) . (46) 


From the numerical calculation we estimate Cw/A ~ 
0.3 — 0.37 in the vicinity of the hrst resonance, which 
can account for the observed reduction of {S) relative to 
5'(G) • Moreover, the shift of Tmax is consistent as well. 

We consider now the second case where several (broad) 
conduction resonances of the barriers are involved, by 
choosing the chemical potential h/Vq = 2.3 and the same 
averaging procedures. The result for (S') can be seen in 
Fig. j^b) which can be compared with the results from 
S^G) in Fig. I^c). Again the characteristic behavior is 
the same in both plots. We see a sign change in the same 
temperature range and also the overall modihcation upon 
increasing a have the same direction. Again the overall 
magnitude is reduced by a factor 3 in (S) with respect 
to S(G). 


B. Fluctuations of (S) 

Within our approach to (S) it is also of interest to ex¬ 
amine the huctuations at low temperature. In Fig. |10| we 
plot the thermopowers corresponding to 50 conhguration 
samples for N = 1000,/r/Vb = 2.3, A = 9. At temper¬ 
atures T/Tq > 10 “^, the thermopower is conhned in a 
very narrow range where the common broad (impurity) 
resonances determine the thermopower’s behavior. As we 
decrease temperature, the thermopowers starts to reflect 
more and more the individual nature of different conhg- 
urations, i.e. the inhuence of the tiny (inter-impurity) 
resonances is increasing. For T/Tq ^ 10“^ the differ¬ 
ent thermopowers are spread over a broad range and 
the averaged thermopower is almost zero, see Fig. ib), 
i.e. the tiny inter-impurity resonances dominate the ther¬ 
mopower’s behavior. 



FIG. 10: A collection of Seebeck coefficients Sq^ as a 
function of temperature which correspond to 50 
different realizations of the underlying parameters for 
N = 1000, ^l/Vo = 2.3, A = 9, a/5 = 0.06. 


In order to characterize the fluctuations in terms of the 
shape of the probability distribution function we com¬ 
puted 50’000 samples of S for 10 “^ < T/Tq < 10“^ 
with N = 2000, ^/Vb = 2.3, A = 9, cr/d = 0 and deter¬ 
mined the mean {S)/So, the standard deviation as/So, 
the skewness {{S — (*S'))^)/cr| and the excess kurtosis 
{{S — {S))'^)/ag — 3 (vanishing for a Gaussian distri¬ 
bution). The results are shown in Fig. |llK a) where we 
skipped the skewness since it fluctuates in a very narrow 
regime (< 0.1) around zero. 

In the temperature range 10“ 2.5 ^ T/Tq <10 the 
small skewness and the vanishing excess kurtosis sug¬ 
gest a Gaussian shaped probability distribution function 
for S. When T/Tq approaches 10“^'^ the standard de¬ 
viation increases and the hnite average {S) decreases 
which demonstrates that the tiny inter-impurity reso¬ 
nances have a growing influence on the thermopower (i.e. 
ksT ^ A) whereas the broad impurity resonances are 
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loosing their impact. 




FIG. 11: (a) Plot of the average {S)/So, the standard 
deviation as/So and the excess kurtosis. (b) 

Comparison of the probability distribution function 
P{S) for T/To = 10“® with a Gaussian and a 
Lorentzian fit. The statistics involves 50’000 samples for 
N = 2000, /i/Vb = 2.3, A = 9, a/6 = 0. 

In the regime where 10“^'^® < T/To < 10“^'^, the 
vanishing average (S) and the finite standard deviation 
reveals that the behavior of the thermopower is entirely 
dictated by inter-impurity resonances {hgT < A). The 
finite excess kurtosis shows that the tails of the probabil¬ 
ity distribution function becomes longer compared to the 
Gaussian distribution function. However, the detailed 
nature of the distribution could not be clarified so far. 

At very low temperatures, 10“® < T/Tq < 10“^-^^, the 
excess kurtosis is strongly enhanced, i.e. the tail must 
have considerably been stretched. In combination with 
the reduction of the standard deviation, we conclude 
that the nature of the probability distribution func¬ 
tion has profoundly changed from a rather Gaussian 
shape to a widely spread shape. In Fig. ©b) we plot¬ 
ted the numerical probability distribution function in¬ 
cluding a Lorentzian and a Gaussian fit on a log-log 


scale. It is evident that the probability distribution is 
in agreement with the Lorentzian distribution function 
for 10“°'® < |S'|/5'o < which is consistent with the 

results by van Langen et aP^. 


V. SUMMARY 

With our study we target some of the basic properties 
of a specially designed disordered finite-size system which 
shows two main features, Anderson localization due to 
random scattering and conductance peaks introduced by 
resonant transmission because of similarly shaped barri¬ 
ers. We show that localization plays an important role, 
as it sharpens the conductance resonances. This in turn 
manifests itself in the enhancement of thermopower. One 
feature is the S oc NT behavior at very low temperatures 
which indicates an enhancement in finite systems with 
stronger localizatiorP^. We also found that the maxi¬ 
mum of I S'I in an intermediate temperature range is again 
proportional to N, if a single resonance dominates the be¬ 
havior. Eventually, the figure of merit could be boosted 
as well by localization effects {ZT oc N), and is only 
truncated by the phonon contribution to the heat trans¬ 
port. 

Within the Landauer-Biittiker formalism we developed 
a procedure based on the harmonic average of the con¬ 
ductance G to determine the mean value of the Seebeck 
coefficient, S(^Qy Examining the validity of this scheme 
in comparison with numerical results, we could show that 
basic results can be trusted, while we could identify the 
short-comings. We also could determine the range of 
validity for Mott’s formula and found a good approxima¬ 
tion for the intermediate temperature range introducing 
a scheme in the expansion in SE, the width of the reso¬ 
nance. 

Although our model has rather severe constraints we 
believe that techniques in producing one-dimensional 
samples with special design might soon reach a stage 
where it could be emulated. We have also demonstrated 
that less perfect scatterers might still give a good perfor¬ 
mance. 
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Appendix A: Computation of the harmonic average 

1. General formula for the harmonic average {G)h 

In this Appendix we address to the problem of calcu¬ 
lating the harmonic average of the conductance {G)h, see 
Eq. (131, for the situation where the transfer matrices of 


the elastic events (scattering and propagation) are known 
and the probability distribution function P decomposes 
into single event probability distribution functions -Pi,p, 
i.e. 




(Al) 


The naive calculation by taking the square of the ab¬ 
solute value of the entity of T and average it afterwards 
normally ends in a recursive formula which might be very 
challenging to solve. Thus, we would like to access the 
problem in a different way and note that 


-Tr 

2 




1 -I- |r|- 


We introduce the harmonic average 


- 1 . 


(A2) 


= J iY[dljddjdVj Pp{lj)P,{6j,Vj) (A3) 


and find by comparing with Eq. (13) that 

2 Go 


{G)h = 


1 -h iTr(ff t) ■ 


(A4) 


The transfer matrix T of a system with N impurities can 
be decomposed such that we can write 




N 


TpilN)Ti{SN, Vn) 




N-1 


fJ{6N,VN)fUlN) . 


(A5) 

We define a mapping / on the vector space of 2x2- 
matrices 


fix) -.= j dl dV dS Ppil)P,i6,V) 

fpmid,V)X f^iS,V)f^il) , 
which allows us to write 

(ff t)^ = / ((ff t)^_i) = ... = /^(l) , 


(A6) 

(A7) 


where we used that (TTl)jv=o = 1- 

The function / is a linear mapping and can therefore 
be fully described by a matrix. We choose the Pauli 
matrices T 2 ,T 3 ,T 4 combined with the identity matrix ti 


Ti = 1 , T2 = 


0 1 
1 0 


TS = 


0 -i 
i 0 


T4 = 


1 0 
0 -1 


(AS) 


as a basis of the vector space. With 

4 4 

A = , Y = fiX) = J2Y,T, 

3=1 3=1 


(A9) 


the linear mapping / is fully described by the 4x4 ma¬ 
trix / = {fij}, 


3 = 1 


(AlO) 


The matrix / is then computed by using Eq. (A6) and is 
given by 


3=1 


(All) 


The matrix fP depends only on the probability distribu¬ 
tion Pp and the transfer matrix Tp. By parameterizing 


rp _ (Til Ti2 
-^P — I rr-i* f 5 


(A12) 


we find that 


r = 




Pp 

0 


0 

0 0 0 

Ip / 


(A13) 


where the 3x3 matrix Fp is 


Pp = J Ppi^) ^ 

( \Tii\^ + \Ti2\^ 2Re(Ti*4ri2) 2 Im(TnTi*2) 

2 Re(TnTi2) Re(r2, + ImiTf, - T^^) 

V-2Im(Tnri2) -MTfi+Tf^) ReiTf^-Tf2), 


(AM) 


and 


lp = I rf/Pp(0(|Tnp-|Ti2p). (A15) 

The matrix /* can be found analogously by applying 
the same scheme on the probability distribution function 
Pi{ 6 , V) and the transfer matrix Tj. 

Since Tr(X) = 2Ai, see Eq. (A9|, and 1 = ri we 
obtain 





(rr) 

= 

[PpP^) 

\ / 

11 



(A16) 


J 11 


Thus, we finally obtain the averaged conductance 

2 Go 


{G)h = 


1 


{hP) 


N 


(A17) 
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2. Computation of {G)h for our model 

We consider a system with the following probability 
distribution functions 


Pi{l) = const 
Py{V)=S{V-Vo) 


Ps{Si) = 


1 




na 


■ exp 


{ 5 , - sy 

2 (t 2 


(A18) 

(A19) 

(A20) 


With the transfer matrices Eq. Q and we find the 
matrices 


1 0 O' 

Fp = I 0 0 0 
,0 0 0 ; 


Fr = 


Q!l + 02 . 


V ••• 


, (A21) 


where 


r r 2 

ai= d6i Ps{5i) cos^(g5i) + — sin^(g(5j) (A22) 
a2 = f dSi Ps{5i) ^sin^(g^i) . (A23) 


We introduce the dimensionless quantities 
E , 2mVod'^ 

and find after some calculations 
sin^ 


(A24) 


ai + a2 = 1 + 


2x(x — 1 ) 


+ q _ c°s(2v'A(x-i)) 

) 4x{x-l) 


Since we have 


{PpF^) 


N 


= {ai + 02 ) 


N 


(A26) 


due to the special shape of Fp, the averaged conductance 
is given by 


{G)h = 


2Gn 


1 + (oi + 02 ) 


N ■ 


(A27) 


Appendix B: The dA-expansion 
1. General formula 

In the limit N —>■ 00, the width of the conductance 
peaks of {G)h goes to zero as 1/^fN. We want to derive 
an expansion in the width of the peak, analogously to 
the Sommerfeld expansion, of the integral 


dx {G)h{x)h{x) , 


(Bl) 


where we use the dimensionless variables x = E 
A = 2mVod^/Ti^ and h is an arbitrary function. 

We concentrate on a single conductance peak at the 
resonance energy Xn = F„/Vb = 1 + (mry/X at perfect 
resonance (cr = 0 ) and expand the averaged conductance 
{G)h, 

{G)h ~ 7 — 7 - — - = Go ■ SnIx - Xn) , 

1 + (1 + ln[x - a;„) 2 ) 

(B2) 

where 


1 


In = 


A3 


(n^)4(l+(^) 


(B3) 


For sufficiently large N, the main contribution to the 
integral Eq. ( |B1[ ) comes from a narrow energy interval 
around Xm which we assume not to be close to 1. Then 
we perform a Taylor expansion of h about Xn and extend 
the lower integral boundary to infinity. The integral is 
then 


dx Sjyix — Xn)h(x) 

dh 


— Aq h{Xn) + Ai 


dx 


A2 d'^h 
2 dx"^ 


with the parameters 

/ OO 

dx 6 n{x - Xn) ■ {x - a;„)\ 

-OO 


(B4) 


(B5) 


The parameter Ai vanishes due to the symmetry of 6^ 
and Ao can be calculated as follows, 


/■°° dx 

(A25) = 


1 


dx 


l+(l + 2#) 


N 


(B 6 ) 


y/N y_oo 1 + exp ( 7 „a; 2 ) A^ 3/2 ’ 

where we applied the identity 

expM = + . (B7) 

The contribution Aq cancels out later and is therefore 
omitted in the following. The integral can be expressed 
in terms of the Riemann zeta function ^ (x) and we obtain 


An 


An 


where 

Ao = -2(y2-l)v^.C(l/2). 

Similarly, we find that 

A 2 


Ao — 


{lnN)^G ’ 


(B 8 ) 


(B9) 


(BIO) 
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where 




(Bll) 


Thus, we get the expansion 


j dx{G)h{x)h{x) « Go ^0 f Kx^) + 


where 


A ^ ^ d^h 
dx'^ 


(B12) 


thermopower in the old notation 

5(G) 

= -5o 


1 




AE A 


, tanh ( —AE ) 

keT 8{kBT)^ U J 


(B16) 


where Sq = ks/e. We approximate 

tanh (/3A_E/2) « sign(£'„ —/i) such that we eventu¬ 
ally end up with the expression 


S(g){T)^-S^ 


(\AE\ A 5E'^ 


V keT 8 {keTf 


•sign(A£;). (B17) 


^ = W =-W77T7W ~ 

Ao 2 v^C(1/2) 


(B13) 


and the full width at half maximum 6x = y/8/jnAf of 
the conductance peak of {G){x). 


2 . Expansion of the thermopower S/^c) 


Then, we apply the expansion on the integrals the 
Cutler-Mott equation Eq. ( [Tl] ) which we rewrite in the di¬ 
mensionless variables x = EJVq, t = T/Tq and /2 = /r/Vb, 


%)=--- 


ksl Jdx {G)h {-dxf) {x - fl) 
e t J dx {G}h {-dxf) 


(B14) 


We find 


S{G) «-T 


ks 1 
e t 


A ^ sr 2 5 , f df 


, (B15) 


where Ax = Xn — p- With r„ = 'Ju/Vq , AE = En — 

Sx = SEjVo and the full width at half maximum 
6E « ^yS/TnN of the conductance peak we obtain the 


3. Expansion of the ZT value 

For a derivation of the electrical conductivity and the 
heat conductivity in the Landauer-Biittiker picture we 
refer to the work of Guttman et aP^. They are given in 
terms of the conductance G by 


^el — AIq , 


(B18) 

(B19) 


where 


K^[G] =JdE G{E) {E - ^i)^ . (B20) 


Then we can apply the 5E-expansion Eq. (B12) and find 
after some calculations that 


Kel = 


AoGoA ,( df 


df 


aei=GoAo{-^j +0{SE^) 


From Eq. (B17) we know that 
ke AE 


S = -- 


e ksT 


0{SE^) 


(B21) 

(B22) 

(B23) 


such that the ZT value for Kph = 0 in the (5i?-expansion 
is given by 


ZT = 
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